SUBSTITUTION PROCESSES IN MOLECULAR EVOLUTION. II. EXCHANGEABLE MODELS FROM POPULATION GENETICS

Substitution processes are of two sorts: origination processes record the times at which nucleotide mutations that ultimately fix in the population first appear, and fixation processes record the times at which they actually fix. Substitution processes may be generated by combining models of population genetics—here the symmetrical-neutral, overdominance, underdominance, TIM, and SAS-CFF models—with the infinite-sites, no-recombination model of the gene. This paper is mainly concerned with a computer simulation study of these substitution processes. The rate of substitution is shown to be remarkably insensitive to the strength of selection for models with strong balancing selection caused by the genealogical drift of mutations through alleles held in the population by selection. The origination process is shown to be more regular than Poisson for the overdominance, TIM, and SAS-CFF models but more clustered for the underdominance model. A class of point processes called Sawyer processes is introduced to help explain these observations as well as the observation that the times between successive originations are nearly uncorrelated. Fixation processes are shown to be more complex than origination processes, with regularly spaced bursts of multiple fixations. An approximation to the fixation process is described. One important conclusion is that protein evolution is not easily reconciled with any of these models without adding perturbations that recur on a time scale that is commensurate with that of molecular evolution.     

Estimation of evolutionary distances between nucleotide sequences

A formal mathematical analysis of the substitution process in nucleotide sequence evolution was done in terms of the Markov process. By using matrix algebra theory, the theoretical foundation of Barry and Hartigan's (Stat. Sci. 2:191–210, 1987) and Lanave et al.'s (J. Mol. Evol. 20:86–93, 1984) methods was provided. Extensive computer simulation was used to compare the accuracy and effectiveness of various methods for estimating the evolutionary distance between two nucleotide sequences. It was shown that the multiparameter methods of Lanave et al.'s (J. Mol. Evol. 20:86–93, 1984), Gojobori et al.'s (J. Mol. Evol. 18:414–422, 1982), and Barry and Hartigan's (Stat. Sci. 2:191–210, 1987) are preferable to others for the purpose of phylogenetic analysis when the sequences are long. However, when sequences are short and the evolutionary distance is large, Tajima and Nei's (Mol. Biol. Evol. 1:269–285, 1984) method is superior to others.     

Contrasting rates of nucleotide substitution in the X-Linked and Y-Linked zinc finger genes

We have sequenced the entire exon (1,180 bp) encoding the zinc finger domain of the X-linked and Y-linked zinc finger genes (ZFX and ZFY, respectively) in the orangutan, the baboon, the squirrel monkey, and the rat; a total of 9,442 by were sequenced. The ratio of the rates of synonymous substitution in the ZFY and ZFX genes is estimated to be 2.1 in primates. This is close to the ratio of 2.3 estimated from primate ZFY and ZFX intron sequences and supports the view that the male-to-female ratio of mutation rate in humans is considerably higher than 1 but not extremely large. The ratio of synonymous substitution rates in ZFY and ZFX is estimated to be 1.3 in the rat lineage but 4.2 in the mouse lineage. The former is close to the estimate (1.4) from introns. The much higher ratio in the mouse lineage (not statistically significant) might have arisen from relaxation of selective constraints. The synonymous divergence between mouse and rat ZFX is considerably lower than that between mouse and rat autosomal genes, agreeing with previous observations and providing some evidence for stronger selective constraints on synonymous changes in X-linked genes than in autosomal genes. At the protein level ZFX has been highly conserved in all placental mammals studied while ZFY has been well conserved in primates and foxes but has evolved rapidly in mice and rats, possibly due to relaxation of functional constraints as a result of the development of X-inactivation of ZFX in rodents. The long persistence of the ZFY-ZFX gene pair in mammals provides some insight into the process of degeneration of Y-linked genes.Correspondence to: W.-H. Li     

Subunits of forming pollen exine and Ubisch bodies as seen in freeze substitutedLedebouria socialis Roth (Hyacinthaceae)

Summary High-pressure freezing/freeze substitution/TEM was employed to investigate anthers of the monocotyledonous angiospermLedebouria socialis Roth (Hyacinthaceae) during early tetrad stage. The initials of the outer sporopollenous pollen wall stratum (=sexine) and of the homologous tapetal products (=Ubisch bodies) are composed of highly regular subunits: clustered globules with a constant diameter of approximately 28 nm. The clusters develop within diffuse accumulations of electron-dense material. This process, interpreted as sporopollenin polymerization, does not necessarily depend on the presence of membrane-bound enzymes. Immunogold labeling with JIM 5 and JIM 7 antibodies revealed that the primexine as well as the dissolving tapetal cell walls, the sites of sexine and Ubisch body formation, respectively, contain un-esterified and methyl-esterified pectins.Abbreviations E-PTA ethanolic phosphotungstic acid - PA periodic acid - UA/Pb uranyl acetate/lead     

Genetic control of amylose content in wheat endosperm starch and differential effects of three Wx genes

The endosperm starch of the wheat grain is composed of amylose and amylopectin. Genetic manipulation of the ratio of amylose to amylopectin or the amylose content could bring about improved texture and quality of wheat flour. The chromosomal locations of genes affecting amylose content were investigated using a monosomic series of Chinese Spring (CS) and a set of Cheyenne (CNN) chromosome substitution lines in the CS genetic background. Trials over three seasons revealed that a decrease in amylose content occurred in monosomic 4A and an increase in monosomic 7B. Allelic variation between CS and CNN was suggested for the genes on chromosomes 4A and 7B. To examine the effects of three Waxy (Wx) genes which encode a granule-bound starch synthase (Wx protein), the Wx proteins from CS monosomics of interest were analyzed using SDS-PAGE. The amount of the Wx protein coded by the Wx-B1 gene on chromosome arm 4AL was reduced in monosomic 4A, and thus accounted for its decreased amylose content. The amounts of two other Wx proteins coded by the Wx-A1 and Wx-D1 genes on chromosome arms 7AS and 7DS, respectively, showed low levels of protein in the monosomics but no effect on amylose content. The effect of chromosome 7B on the level of amylose suggested the presence of a regulator gene which suppresses the activities of the Wx genes.     

Role of the Cys 2-Cys 10 disulfide bond for the structure, stability, and folding kinetics of ribonuclease T1.

The Cys 2-Cys 10 disulfide bond in ribonuclease T1 was broken by substituting Cys 2 and Cys 10 by Ser and Asn, respectively, as present in ribonuclease F1. This C2S/C10N variant resembles the wild-type protein in structure and in catalytic activity. Minor structural changes were observed by 2-dimensional NMR in the local environment of the substituted amino acids only. The thermodynamic stability of ribonuclease T1 is strongly reduced by breaking the Cys 2-Cys 10 bond, and the free energy of denaturation is decreased by about 10 kJ/mol. The folding mechanism is not affected, and the trans to cis isomerizations of Pro 39 and Pro 55 are still the rate-limiting steps of the folding process. The differences in the time courses of unfolding and refolding are correlated with the decrease in stability: the folding kinetics of the wild-type protein and the C2S/C10N variant become indistinguishable when they are compared under conditions of identical stability. Apparently, the Cys 2-Cys 10 disulfide bond is important for the stability but not for the folding mechanism of ribonuclease T1. The breaking of this bond has the same effect on stability and folding kinetics as adding 1 M guanidinium chloride to the wild-type protein.     

A scanning electron microscope study of normal and vitrified leaves from Datura insignis plantlets cultured in vitro

The surface anatomy of normal and vitreous leaves of plantlets obtained from Datura insignis Barb Rodr nodal segment cultures was compared using scanning electron microscopy. Normal and vitrified leaves are similar in several ways. They are both amphistomatic, and have similar distributions of glandular and non-glandular trichomes. Stomata have similar length, diameter and distribution in normal and vitreous plants. Immature stomata, which have closed pores, and plugged stomata, which contain an amorphous material between their guard cells, occur in both normal and vitrified leaves. Normal and vitreous leaves differ in the frequency of normal and abnormal stomata. Normal stomata have kidney-shaped guard cells and resemble closely those found in field-grown plants, whereas abnormal stomata have deformed guard cells. Normal stomata represent approximately 80% of the total number of stomata in normal leaves, but only 7% of the total number of stomata in vitreous leaves. Abnormal stomata represent 90% of the total number in vitreous leaves. The deformation of guard cells could possibly be a mechanical impediment to stomatal function.     

载脂蛋白多基因家族分子进化的研究

与脂质运输有关的载脂蛋白基因构成一个复杂的多基因家族。为探讨这种演化时间长的基因家族的进化规律,本文首先建立了一种在非均衡进化速率条件下计算系统发生树中任意分支长度的简易方法,并可在此基础上算出无根分支系统树中分歧年代的期望值。进一步对本文科１０个种属共２６种载脂蛋白的系统演作作了实际分析,结果提示：①ＡｐｏＡ－Ｉ＇ＡｐｏＡ－ＩＶ,ＡｐｏＥ及ＡｐｏＡ－ＩＩ的共同祖先可能在奥陶纪水生脊椎动物中就已存